Possible new phenomenology in turbostratic graphene bilayer 
as a solid-state realization of a two-brane universe model 
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It is suggested that two turbostratic graphene layers can be described through a formalism previ- 
ously introduced for the study of braneworlds in high energy physics. It is demonstrated that any 
graphene bilayer can be conveniently described as a "noncommutative" two-sheeted (2+1) spacetime. 
The model Hamiltonian contains a coupling term connecting the two layers which is a reminiscent of 
the coupling which may exist between two braneworlds at a quantum level. In the present case, this 
term is related to a K — K' intervalley coupling. Phenomenological consequences are emphasized 
such as exciton swapping between the two graphene sheets. An experimental device is suggested to 
study this possible new effect. If demonstrated, this effect could be an indirect way to test some 
theoretical concepts of the braneworld hypothesis. Technological applications are underlined. 

PACS numbers: 72.80.Vp, 78.67. Wj, 11.25. Wx, 02.40.Gh 



I. INTRODUCTION 



During the last few years, graphene has taken a grow- 
ing importance in solid states physics [H-I22]] . Indeed, it is 
an amazing case of two-dimensional carbon crystal, and 
its remarkable properties make it a strategic material for 
future nanotechnologies. From a fundamental and theo- 
retical point of view, graphene is also a solid-state real- 
ization of a (2 + 1) spacetime in which massless fermions 
live. Thus, graphene is a remarkable tool to experiment 
some theoretical concepts of low-dimensional electrody- 
namics and quantum dynamics flll - ilol ]. Graphene could 
be also an efficient tool to test some theoretical ideas tak- 
ing their origins in high energy physics. During the last 
two decades, the possibility that our observable (3 + 1)- 
dimensional Universe could be a sheet (a 3— brane or 
braneworld) embedded in a (N + l)-dimensional space- 
time (called the bulk, with N > 3) has received a lot 
of attention (23|. Such an exotic concept appears very 
productive to solve puzzling problems beyond the stan- 
dard model of particles (23|. Nevertheless, braneworld 
models remain complex especially to confront theoretical 
solutions with experiments. However, in recent papers, 
a mathematical and physical tool was introduced to sim- 
plify the study of a universe made of two branes (2^ - |27j . 
At low energies, the quantum dynamics of Dirac fermions 
in such a two- world can be rigorously described by a more 
simple and equivalent model that corresponds to a two- 
sheeted spacetime in the formalism of the noncommuta- 
tive geometry [24l HH. The model reveals a coupling term 
in the Hamiltonian which connects the two branes: Par- 
ticles endowed with a magnetic moment undergo Rabi 
oscillations between those worlds under the influence of 



'Electronic address: michael.sarrazin@fundp.ac.be 
t Electronic address: |f.petit6ibcrc.be| 



a magnetic vector potential [244271 . 

On another hand, recent works on graphene underline 
the importance of electronic transport in turbostratic bi- 
layers [3|-[9j ■ Since a graphene sheet can be considered as 
2— brane embedded in a (3 + 1) bulk, a graphene bilayer 
could be a solid-state realization of a universe contain- 
ing two branes (a two-brane universe). In the present 
paper, we show that this analogy is well-sounded and 
demonstrate the possibility to apply our noncommuta- 
tive geometry tool to study a graphene bilayer. This 
approach suggests that exciton swapping may occur be- 
tween the two layers. Thus a possible new experimental 
phenomenon would be a solid-state counterpart of the 
particle oscillations predicted in brane theory [23 - |27| . 

In section UH we recall the basic assumptions under- 
lying the description of electron and hole in graphene 
through a Dirac equation formalism. Next, in section 
THl we present our model of fermion dynamics in a two- 
sheeted spacetime and its adaptation to describe a set 
of two graphene layers. Then, in section IIV1 using a 
tight-binding approach, we show that two turbostratic 
graphene layers are a prerequisite to get a K — K' inter- 
valley coupling between two perfect graphene layers in 
mutual interaction. This coupling may result in exciton 
swapping between the two layers as shown in section fVl 
Finally, in section lVll an experimental device is suggested 
to investigate this new effect. 



II. GRAPHENE ELECTRONIC PROPERTIES 



Graphene is a one-atom thick layer made of sp 2 car- 
bon atoms in an hexagonal lattice arrangement (Fig. la) 
[3, Q . Self-supported ideal graphene is a zero-gap semi- 
conductor. In the vicinity of the six corners (called Dirac 
points) of the two-dimensional hexagonal Brillouin zone 
(Fig. lb), the electronic dispersion relation is linear for 
low energies (Fig.lc). Electrons (and holes) can then 
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FIG. 1: (a) Hexagonal lattice of graphene with the two sub- 
lattices A and B. ai and ai are the vectors of the unit cell, (b) 
Brillouin zone of the hexagonal lattice, (c) Energy behavior 
in the vicinity of the Dirac points K and K' . 



It can be noticed that the above (2 + 1)-Dirac equation 
can be easily extended to its (3 + l)-dimensional version. 
7 3 and 7 5 matrices (such as 7 5 = i7°7 1 7 2 7 3 ) can be 
introduced and we may consider for instance: 
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The Clifford algebra is verified since {7^,7"} = 
{l 5 ,l"} = and H7 5 ) 2 = -1, where g» v is the 
four-dimensional metric tensor of the Minkowski space- 
time (/i, v = 0, 1, 2, 3). Note that the j 3 and 7 s matrices 
are interchangeable through substitutions 7 3 — > ij 5 and 
— i-f 5 — > 7 3 which lead to equivalent descriptions. More- 
over, it is well known that 7 s can be also used to define a 
five-dimensional Dirac equation as shown in section [nil 



be described by a Dirac equation for massless spin— 1/2 
particles in an effective (2+1) spacetime [2f . While mass- 
less Dirac fermions propagate at the speed of light in the 
(3 + 1) Minkowski spacetime, in graphene the effective 
massless Dirac fermions propagates at the Fermi velocity 
(vp w 10 6 m-s" 1 in the present case). On a graphene 
layer, the Hamiltonian of the effective Dirac equation is 
given by Q: 



H± = —ifiVF((j\d x ± (Tidy) 



>j f a 3 



(1) 



where "+" (respectively "— ") refers to the K (respec- 
tively K') point of the Brillouin zone of the graphene 
hexagonal structure (Fig. la), (k — 1,2,3) are the 
usual Pauli matrices. For a self-supported graphene sheet 
the mass term m is equal to zero and electrons (and 
holes) behave as relativistic quasiparticles. Neverthe- 
less m may differ from zero in the case of a sheet de- 
posited on a substrate 16]. Using m — > mvp/fr and 
(xq, x\, x 2 ) = (vFt,x,y), from Eq. (QJ it is possible 
to conveniently describe the electron (or hole) dynam- 
ics through an effective Dirac equation such that |2j : 



{i'f ri d v — m)tp = 



(2) 



with 77 = 0, 1, 2 and 



o / (73 \ 1 / ia 2 2 / -ia t 
7 ^ a 3 ) ' 7 ^ ia 2 ) ' 7 { ia x 

(3) 

such that {7 , ',7 1? } = 2gi*> (77,1? = 0,1,2) with 
gird =diag(l, — 1, — 1). The wave function is defined as 

ijj = ( ^ ^ where x (respectively 0) is related to the 
wave function on K (respectively K'). In addition, \ 
(respectively 9) can be written as x = f ^ 1 (respec- 
tively 9 = ( ^ A ^ ) where A and B are related to the two 
sublattices of the graphene sheet (see fig. la). 



III. TWO-LAYER GRAPHENE AS A 
"NONCOMMUTATIVE" TWO-SHEETED 
SPACETIME 



Let us consider a graphene layer as a 3-brane, i.e. a 
three-dimensional space sheet, for which one dimension 
(say x 3 ) is reduced to zero. We suggest to derive the 
graphene bilayer system description from the two-sheeted 
spacetime model introduced in previous works [23 - [27j by 
making X3 — > 0. 

In a prior work, the relevance of the two-sheeted ap- 
proach was rigorously demonstrated for braneworlds de- 
scribed by domain walls Indeed, when one studies 
the low-energy dynamics of a spin— 1/2 particle in a two- 
brane Universe, the quantum dynamics of this particle is 
equivalent to the behavior it would have in a two-sheeted 
spacetime described by noncommutative geometry [24j |. 

Specifically, a two-sheeted spacetime corresponds to 
the product of a four-dimensional continuous manifold 
with a discrete two-point space and can be seen as a five- 
dimensional universe with a fifth dimension reduced to 
two points with coordinates ±<5/2. Both sheets are sep- 
arated by a phenomenological distance 5, which is not 
the real distance between the graphene layers as shown 
in the next section. Mathematically, the model relies on 
a bi-euclidean space X = M4 x Z 2 in which any smooth 
function belongs to the algebra A = C°°(M) © C°°(M) 
and can be adequately represented by a 2 x 2 diagonal 
matrix F =diag(/i, f 2 ). In the noncommutative geom- 
etry formalism, the expression of the exterior derivative 
D = d + Q, where d acts on M4 and Q on the Z 2 internal 
variable, has been given by Connes [23]: D : (fi,f 2 ) — ► 
(df 1 ,df 2 ,g{f 2 - /i),.g(/i - f 2 )) with g = 1/6. Viet and 
Wali [29( have proposed a representation of D acting as a 
derivative operator and fulfilling the above requirements. 
Due to the specific geometrical structure of the bulk, this 
operator is given by: 
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H = 0,1,2,3 and D 5 
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(5) 

where the term g acts as a finite difference operator along 
the discrete dimension. Using (5), one can build the 
Dirac operator defined as = T N D N = T»D U + T 5 D 5 . 
It is then convenient to consider the following extension 
of the gamma matrices (by using the Hilbert space of 
spinors [28}): 
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(6) 



In the present work, 7^ and 7° = ij""f'"fj" 
Dirac matrices defined by relations © and ^ relevant 
for graphene. We can therefore introduce a mass term 
M = ml_8x8 a s in the standard Dirac eq uation. The 
two-sheeted Dirac equation then writes |24l - t26| : 



(0 - M) W = (iT 



N 
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(7) 



m 



= 



with 



i>0 



the two-sheeted wave function. In this 



notation, the indices "a" and "/3" discriminate each sheet 
[2^ - |26j . i.e. each graphene layer when x% — ► 0. Each 
component of the wave function tp is then the probabil- 
ity amplitude of the electron (or hole) in each graphene 
sheet. It is important to point out the Lagrangian term: 



(8) 



which ensures the coupling between each graphene layer 
through K~K' processes as explained in section llVl This 
C c term is the main reason for this paper as it allows 
particle swapping between the graphene layers. 

Let us now introduce the effect of an electromag- 
netic field, i.e. an U{1) gauge field. To be consis- 
tent with the two-sheeted structure of the Dirac field 
'F in Eq. 0, the usual U(l) electromagnetic gauge 
field should be replaced by an extended U(l) ® U(l) 
gauge field [HrHl]- Nevertheless, in the present work, 
we assume that electromagnetic field sources are out of 
the graphene layers. The group representation G = 
diag(exp(— iqA a ), exp(— iqAp)) is therefore reduced to 
G = diag(exp(—iqA),exp(—iqA)). We are looking for 
an appropriate gauge field such that the covariant deriva- 
tive becomes Jj) — > with the gauge transformation 
rule 4- = GA'G^ — iG [Jfi,u rac , G t ] . A convenient choice 
is (H-H 



tqj 







(9) 



(respectively A@) is the magnetic vector potential A^ 
on the graphene layer a (respectively /3). According to 
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FIG. 2: (Color online) Sketch of the two turbostratic graphene 
layers under consideration. Both sheets are rotated with re- 
spect to each other with an angle 6 ~ 21.787°. t\ and £2 are 
the vectors of the Moire unit cell. 



the appropriate covariant derivative, the introduction of 
the gauge field in Eq. © leads to [24j2r| 



i9T 



*57 5 \(^a. 
i^(d u + iqAP)-m J \ 

Of course, for graphene sheets, we have X3 — which 
corresponds to two bidimensional sheets instead of three- 
dimensional space sheets. In addition, we will assume 
that the experiment under consideration hereafter (sec- 
tion lVI[) is such that A^ is parallel to each graphene sheet 
(A 3 = 0). 



IV. K — K' COUPLINGS IN TURBOSTRATIC 
GRAPHENE LAYERS 



Let us now suggest a heuristic argument to justify the 
use of our brane model for studying the graphene proper- 
ties. A full demonstration of this statement is likely well 
beyond the scope of this paper but simple arguments suf- 
fice to be convinced by the analogy graphene bilayer/two- 
sheeted spacetime. In braneworld models, we have sim- 
ply to consider the interaction between one fermion and 
domain walls described by a scalar field 24[. By contrast, 
a bilayer graphene is formally a many-body problem. 
Therefore we should normally consider the whole dynam- 
ics of carbon atoms and their electrons. This would be a 
very complicated task of course. Anyway, the existence 
of coupling terms proportional to g is straightforward for 
turbostratic graphene layers. When two graphene layers 
are twisted with respect to each other, a typical Moire 
pattern can be observed (Fig. 2). This occurs when 
both layers are commensurate, i.e. when two specific 
kind of atoms of each layer can be superimposed peri- 
odically The Moire pattern can be then described 
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through a periodic unit cell defined by vectors ti and t 2 
(see Fig. 2) and can only exist for a specific rotation angle 
9 = 6 Pt g (with p, q G PJ) between both layers. 

Let us define ai = <2o(l/2, y3/2) and a 2 = 
ao(— 1/2, V3/2) the vectors of the real space which define 
the unit cell of the first graphene layer (see Fig. la), do is 
the lattice parameter. Two kinds of commensurate struc- 
tures can be considered [3, 0] • The first one is such that 
the vectors of the Moire unit cell are t x = p&i + (p + q)a 2 
and t 2 = — (p + q)ai + (2p + q)& 2 such that gcd(q, 3) = 1. 
The second case is such that ti = (p + q/3)ai + (q/3)a 2 
and t 2 = — (q/3)ai + (p + 2q/3)a 2 with gcd(q, 3) = 3. In 
both case, the rotation angle p ^ q between both sheets is 
given by 0,11]: 



cos 9 p ^ = 



3p 2 + 3pq + q 2 /2 
3p 2 + 3pq + q 2 



(11) 



In the first layer, the first K Dirac cone is located at K = 
(47r/(3a ))(l, 0) while the K' Dirac cone is at K' = -K. 
By contrast, in the second layer, due to the rotation the 
K Dirac cone is located at K e = (47r/(3ao))(cos sin#) 
whenever the K' Dirac cone is at K 6 ' = -K e [!, [|. 
Let Gi and G 2 be the vectors of the unit cell of the 
reciprocal lattice of the Moire pattern. Obviously, the 
Moire pattern can be responsible for coupling between 
valleys of each layer [l-Q- Indeed, we get 



G = K-K e = -(K'-K' e ) 
for K — K couplings, and 

G C = K-K' 9 = -(K'-K e ) 



(12) 



(13) 



for K — K' couplings. When gcd(g, 3) = 1, then G = 
-(?/3)(2Gi + G 3 ) and G c = -(2p + q)G 2 . While, 
when gcd(g, 3) = 3, then G = — (q/3) (Gi + G 2 ) and 
G c = (l/3)(2p + q) (Gi - G 2 ) . The greater G and G c 
are, the weaker the couplings are. As a consequence, 
one should consider the lowest values of p and q. A 
similar consideration leads us to expect that K — K in- 
terlayer couplings are usually stronger than the K — K' 
ones. Then, for the purposes of our study, it should be 
relevant to consider a structure which can suppress the 
K — K couplings while enhancing the K — K' interlayer 
couplings. We may consider for instance the case such 
that gcd(q, 3) = 1 with q = 1. Indeed, in that case 
G = —(1/3) (2Gi + G 2 ) is not a vector of the recipro- 
cal lattice. By contrast G c = — (2p + 1)G 2 is always a 
vector of the reciprocal lattice and is such that G c ps 2K 
whatever p. The first relevant value to be considered is 
then p = 1. In this case, 6^1 s» 21.787° and we obtain 
the specific structure shown in Fig. 2. 

In a tight-binding approach it is possible to define the 
operator d> a ^ (respectively a a (p),j) which creates an 
electron (respectively a hole) on the site j of the sublat- 
tice "A" on the a graphene layer (or on the /3 graphene 
layer). The same convention is used for the sublattice 



If one considers the interlayer coupling, one gets 
for the turbostratic system 



where the energies t uv j (with u — A, B and v = A, B) 
are related to the interlayer hopping between the nearest 
sites of each layer. This dependence of t uv j vs. the 
location j is very specific for two turbostratic graphene 
layers. In the structure considered here, we can see that 
no AA site exists by contrast to the AB sites (Fig. 2). 
We then assume that tAA,j ~ ^BB,j ~ 0. In addition, 
t AB (Rj) = t AB ,j = -H when Rj = (2/3)(ti + t 2 ) + 
(nti + mt 2 ) (with ti = ai + 2a 2 and t 2 = — 2ai + 3a 2 ) 



and tsA^B-j) = ^ba, 



-t 1 when K 3 = (1/3) (tj 



(nti +mt 2 ), with n, m £ N. tABj and tsAj are equal to 
zero elsewhere. We use the following Fourier transform 
of the operators: 



a a(0){ r j) - a a(l3),j 



E 

k 



1 



N 



»(/3),qj 



(15) 



with a similar convention for b a tp\j and t uv _j 



are 

vectors of the reciprocal lattice of graphene sheet a (or 
/3). Let us consider a single particle state with momen- 
tum k such that we can consider the restricted Fourier 
representation of the Hamiltonian: H c = i/ Cj K+k + 
H c ,K'+\t + ff c .K«+k + -ff c ,K e '+k snch that H c — ^% c ^ 
with 



He = -t'a x 



and 



= ( a a .K h a .K a a ,K' b a> K' 

~ ( i>i V>| ) 



(16) 



,K bp : K a/3,K> bfj : K' ) 

(17) 



Using the above notations, and after a convenient SU(2) 
rotation such that b ib, the Lagrangian term related to 
T~L C in Dirac notation becomes C c = ^iT b D^ , i.e. Eq. 
© related to Eq. ©. 

Now, the coupling constant g can be then defined as 
g = t' /hvp and the phenomenological distance is S = 
HvF/t'. Basically, g and S must depend on the real dis- 
tance d between each graphene sheet. Indeed, the hoop- 
ing energy t' varies as Jp: t' ~ to exp(5.43 • (1 — d/a m )), 
with t ~ 0.3 eV 0, U HI, and here d is the distance 
between two layers, while a m is the nearest interlayer 
distance, a m = 3.35 A. For closest layers (d = a m ), we 
get S of about 22 A (i.e. g ~ 4.5- 10 s m^ 1 ). As an indica- 
tion, note that for d — 2et m (respectively d = 5a m ), one 



gets g w 2 • 10 m (respectively g « 1.7- 10 
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V. PHENOMENOLOGY OF THE MODEL 



Following previous works [24l426| , we focus on the non- 
relativistic limit of our Dirac like equation. Defining 
V = (di,d 2 ), A = (A 1 ,A 2 ), a = (a u a 2 ) and B 3 = 

d x A 2 - d 2 A x and using: F A(B) = ( X a A(B) ) ,and fol- 

V °MB) J 

lowing the well-known standard procedure, a two-layer 
Pauli equation can be derived from Eq. (|T0|) [23 - [27j : 

4fe) = < H °+ H ™>fe) (18) 

where Fa , a and Fa,p correspond to the wave func- 
tions in the graphene layers a and (3 respectively. The 
Hamiltonian Ho is a block-diagonal matrix such that 
Ho =diag(H Q , H^), where each block is simply the effec- 
tive Pauli Hamiltonian expressed in each graphene layer 

ft?" / Q \ ^ 

H a (/3) = [y - *^ A q(/3)J +M3#3,a(/3)+ 1 4(,3) (19) 

such that A Q and Ap correspond to the magnetic vector 
potentials on the layers a and (3 respectively. The same 
convention is applied to the magnetic fields B Q (^) and to 
the potentials V a (p) . In the following, since we consider 
neutral excitons, we can set V a (g\ = 0. In addition, we 
will show hereafter that B 3 a ^ = in the device under 
consideration (see section |Vl|). We set /i = j(ft/2)a where 
7 is the pseudo-gyromagnetic ratio and fi the pseudo- 
magnetic moment related to the pseudo-spin of the par- 
ticle [23] . With this choice, the present approach can be 
extended to any particle endowed with a magnetic mo- 
ment whatever its spin value. 

In addition to these usual terms, the two-layer 
graphe ne Hamiltonian comprises also a new specific term 

^ igii ■ {A a — A^} ) y ' 

H cm is obviously not conventional and describes the cou- 
pling of the layers through electromagnetic fields. It van- 
ishes for null magnetic vector potentials. Intuitively, the 
coupling generated by this term will imply Rabi oscilla- 
tions of electrons or holes between both graphene sheets 
due to electronic derealization. 



VI. EXCITON SWAPPING BETWEEN TWO 
GRAPHENE LAYERS AND EXPERIMENTAL 
DEVICE 



Guided by the previous equations, we now suggest an 
experimental approach for testing exciton swapping be- 
tween two graphene layers. An incident electromagnetic 




FIG. 3: (Color online). Sketch of a feasible experimental 
setup, (a): Basic setup. Two coaxial annular magnets, with 
different inner and outer diameters, covered by an insulating 
material. The upper ring is filled up with an opaque material. 
Two magnetic fields (B a and Bp) turn around the symmetry 
axis of the magnets and are opposite. Two graphene layers 
(G a and Gp) are considered, each one deposited on a face of 
a magnet. The geometry of the device allows the existence of 
two opposite magnetic vector potentials (A a and Ap), each 
one in the vicinity of a graphene layer. An incident photon 7i 
pumps an exciton e on G a . A photon j-t resulting from the 
exciton decay on Gp can be recorded, (b): Full setup. Rect- 
angular array of annular devices deposited on a transparent 
substrate (blue layer). The area between the toroidal mag- 
nets is filled with an opaque material (yellowish layer). Such 
a setup allows to enhance the recorded signal by increasing 
the graphene area. 



wave with a convenient energy can excite an electron-hole 
bound pair (i.e. an exciton) [l7l - [2l| on a first graphene 
layer (G a ). In the best of our knowledge, studies re- 
lated to the magnetic moment of exciton in graphene are 
still lacking. Nevertheless, exciton should exhibit reso- 
nance states endowed with non-zero magnetic moment fi 
[30| due to the combination of the electron/hole magnetic 
moments [22| , possibly supplemented by an orbital mag- 
netic moment. One can then expect to induce a coupling 
through Hem between G a and a second graphene layer 
Gp leading to a swapping of the exciton from G a towards 
Gp. Afterwards, the exciton decay on the second layer 
could be recorded. 

The required magnetic vector potentials can be pro- 
duced with the following device. Let us consider two 
coaxial annular magnets covered by an insulating ma- 
terial (see Fig. 3a). Both magnets have the same rect- 
angular section. Both magnetic fields (B„ and B^) in- 
side the magnets turn around the symmetry axis of the 
magnets and are opposite. B a = Bp = outside the 
magnets due to the toroidal topology [3l|. Only a mag- 
netic vector potential A exists outside the magnet [31( 
(i.e. V x A = 0). Boundary conditions result from 
§ c A • d\ =<f>, where C is a contour on a magnet (see 



6 



Fig. 3) and $ the magnetic flux inside a magnet. The 
geometry of the device leads to two opposite magnetic 
vector potentials (A Q and Ap), each one in the vicinity 
of a graphene layer (G a and Gp) deposited on a face of a 
magnet (see Fig. 3a). A straightforward calculation shows 
that |A Q — A/3 1 ~ 2Aod/ (L + I), with d the distance be- 
tween the two layers, L and / are the length and width 
of rectangular section of the magnets. If one considers 
a superconducting magnet, then A ~ nh / (Ae(L + I)) , 
where n is an integer (h is the Planck constant and e the 
electric charge), due to the magnetic flux quantization 
[32| . For instance, if L = 1 pm and I — 10 nm [3l[ and 
with d = 2a m , one gets |A Q — Ap\ ss 1.4 • 10~ 12 T-m for 
n = 1. 

The insulating material is the substrate on which the 
graphene layers are deposited. This allows a gated 
graphene leading to electron and hole that share the same 
effective mass [161 ] . The efficient graphene area can be in- 
creased by using a large array of micro-annular devices 
(see Fig.3b). 

The excitonic swapping can be described as follows. 
One looks for an exciton wave function in the form: 



!*(*)} = 



F A , a {t) 
Fa A*) 



(21) 



a a (t) 







a p{t) 







where it is assumed that fi^ s = ±^i$ s , i.e. fy s is an 
eigenstate of with an eigenvalue p different from zero. 
For an exciton, the lowest expected value can be esti- 
mated by p ~ eh/m [13], i.e. p « 3.5 • 10" 22 J-T" 1 for 
an effective electron/hole mass about 0.3 eV [l6|]. Putting 
Eq. (|21[) into the Pauli equation (fT8"]) leads to the follow- 
ing system of coupled differential equations: 



—a a = -Ka -{l/2)T o a a + 5(t-ti) (22) 



and 



— ap = na a - (l/2)F ap 



(23) 



with k = pg (A a — Ap) /h. With the above mentioned 
values, one can roughly estimate k « 2.1- 10 9 rad-s -1 . To 
is the exciton decay rate conveniently introduced in the 
equations in agreement with the lifetime r of the exciton 
(ro = t _1 ). We assume that r is included between 10 fs 
and 200 ps HJHj (5-10 9 s" 1 < r < 10 14 s" 1 ). 6(t-U) 
is a Dirac delta source such that the exciton is created 
at t = ti in the layer a. Then, a a (t = ti) = 1 and 
ap{t = ti) = 0. The number of excitons is then given 
by M a = J2i a *a a a (respectively Np = J2i °*p a fi) in la y er 
a (respectively in layer /3). In the continuous limit such 
that M. excitons are produced per second, from Eqs. (|22|) 
and (|2"3"j) . one easily obtains three Bloch-like equations: 



dt 



-kU - T a Af a + M 



(24) 



and 



and 



■Mb = kU — T n Mp 



dt 



—U = 2nM a - 2nMp - T U 



(25) 



(26) 



with U = J2i( a a a P + a aO*p) ■ Since layer a is continuously 
supplied with new excitons thanks to an incident photon 
flux To, the exciton source is such that M. = p c ff2<). p e e is 
the photon-to-exciton conversion efficiency. Eqs. (|24|) to 
must present short-time transient solutions due to 



— ToM a (p-) and — TqU terms. As a consequence, we look 
for stationary solutions such that dM a /dt — dMp/dt — 
dU/dt = 0. Eqs. (24]) to (26]) can be then trivially solved. 
The number of excitons in each graphene layers are: 



M a 



2k 2 



T (4« 2 + r 2 ) 



M, and Mb 



2k 2 



r (4« 2 + r 2 ) 



M 

(27) 

and the number of new created excitons balances the 
number of decaying excitons, i.e. M = r (M a +Mp). 
Note that in the present approach, we do not consider 
any saturation effect regarding to the number of excitons 
per unit area. Then for a fixed area, M a + Mp should 
be limited and MI/Tq likewise. As a consequence, for a 
given value of Ml , the present approach is not valid when 

r -> o. 

The photon flux I t emitted from the second graphene 
layer /3 is It = nT^Mp where n is the number of photons 
that results from the exciton decay. The effective optical 
transmission coefficient T of the device is T = Tt/Io, and 
one gets: 



T = np cS 



2k 1 



4^ 



r 2 , 



(28) 



The excitons transferred from layer a to layer /3 are then 
detected through recorded photons due to excitonic decay 
(see Fig. 3a). Let us consider the simplest process such 
that np c fi = 1, i.e. every exciton decays into a single 
photon, and each photon creates a single exciton [34] . 
With the above values, the best expected transmission T 
could reach 21 %, which is of course a fair value in an 
experimental context. 



VII. CONCLUSIONS 

Using a theoretical approach previously considered to 
describe a Universe made of two braneworlds [UHS] , we 
have proposed a new theoretical description of the phe- 
nomenology of two turbostratic graphene sheets. This 
suggests a new way to describe multilayer graphene, 
which will deserve further studies. We have shown that 
the transfer of excitons between the two graphene sheets 
is allowed for some specific electromagnetic conditions. 
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While the excitons are produced by incident light on the 
first graphene layer, photons could be recorded in front 
of the second graphene layer where the swapped exciton 
decays. Such a system is a solid-state realization of a 
two-brane Universe, for which it has been shown that 
matter swapping between two braneworlds could occur 
[2^ - |27l |. As a consequence, any experimental evidence 
of this effect in graphene bilayers would be relevant in 
the outlook of braneworld studies. Since the magnetic 
fields inside the magnets can be controlled with a tran- 
sient external magnetic field, it is possible to turn on or 



off the device. We can then expect to get a new kind of 
electro-optic light modulator with hysteresis. 
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